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Abstract 

In this paper, an efficient algorithm of logarithmic transformation to Hirota bilinear form of the 
KdV-type bilinear equation is established. In the algorithm, some properties of Hirota operator 
and logarithmic transformation are successfully applied, which helps to prove that the linear 
terms of the nonlinear partial differential equation play a crucial role in finding the Hirota bilinear 
form. Experimented with various integro-differential equations, our algorithm is proven to be 
more efficient than the algorithm referred by Zhou, Fu, and Li in getting the Hirota bilinear form, 
especially in achieving the coefficient of the logarithmic transformation. 

Keywords: Hirota bilinear form, KdV-type bilinear equation, D-operator, Logarithmic 
transformation. Symbolic computation. 



1. Introduction 



In 1971, Hirota developed a direct method, Hirota bilinear method, to construct the exact 
multi-soliton solution of integrable nonlinear partial differential equation (NPDE) [1-3]. Once 
the Hirota bilinear form(HBF) of a NPDE is given, there are several ways to solve it. The solu- 
tions of NPDE can be constructed from HBF by the typical perturbation expansion method j!]] 
and others [4,5]. Therefore, the key step of the Hirota direct method is to transform the NPDE 
into the Hirota bilinear form (HBF). Hietarinta designed a program on searching for integrable bi- 
linear equations such as the KdV-type [Ql, mKdV-type SG-type Ht], NLS-type |@] equations. 
In 1992, Hereman and Zhuang jlOll gave a summary on types of bilinear equations. In recent 
years, some algorithms for generating bilinear form of NPDE are described and Maple packages 



Bilinearization and HBFTrans are established by Zhou, Fu, Li 111 IL 11211 and Yang, Ruan II13I1 . 
respectively. The package Bilinearization can construct the HBF of many NPDEs by solving 
a system of over-determined algebraic equations with respect to the combinatorial coefficients. 
However, the general ansatz of the bilinear form of NPDE in Ref. ll lUl2fl . which relies on the 



WTC method Ill4tl and the HBM method lllSll . is so complicated that the efficiency of perfor- 
mance is relatively low. Without depending on these two methods, our method can obtain the 
HBF of the KdV-type equations within shorter time. In this paper, some properties of Hirota- 
operator are taken advantage of, which brings forward to a more efficient algorithm for finding 
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the HBF of KdV-type equations in NPDE. We take a series of classic KdV-type equations for 
instance, to demonstrate the validity of our algorithm. Furthermore, the implementation of the 
algorithm in Maple is applied to automate the tedious computation for the construction of the 
HBF of KdV-type equations in NPDE. 

2. Hirota Bilinear Method 

Now, let us briefly review the Hirota derivatives. In 1971, Hirota developed the Hirota deriva- 
tive, which is also called D-operator iJ]. For 1 + 1 dimensions, the D-operator is defined by 

■ 8 = - dA\d, - dr)'"f{x, t)g{x', m, « = 0, 1, 2, 3, • ■ ■ . (2.1) 

where /(x, t) and g{x, t) are dififerentiable functions of x and t, respectively. 
From the definition, there are some properties ^ of the D-operator, 

DlD'^f ■ g = D"/D"J ■ g = D"-'DTDJ ■ g (2.2) 
D"^D';'f ■ 1 = d",d';'f, where d", = d"/dx". (2.3) 

D:.Drf-g^{-iy"^"D:D';'g-f (2.4) 

DlD'I'f ■ / = 0, if m + n is odd; (2.5) 



furthermore, for arbitrary independent variable x, a recursive definition of D-operator is. 



I fR, for yi -Q 

D"J-g^Y^\ , (2.6) 

\D"-'f,-g-D"-'f-g,, for «>0. 



As an example, the Korteweg de Viries (KdV) equation 

Uxxx + 6mm.v + m, = (2.7) 
where u - u{x, f), can be transformed through the dependent variable transformation, 

M = 2(ln/),, (2.8) 

into 

<^ I f.af - fx ft + fxxxxf — ^f XXX fx + 3/„ i _ q (2 9) 



dx\ /2 

from which the bilinear equation below is obtained 

/«/ - + fxxxxf - 4fxxxfx + 3fl = Cf\ (2.10) 

where c is a constant of integration. Equation ( 12.10b . with c = 0, may also be written concisely 
in terms of D-operators as 

(Dt + DxD,)f-f^O. (2.11) 

Thus, the Hirota bilinear form of the KdV equation is obtained. 
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3. Principles Of Our Algorithm 

3.1. The Properties ofHirota Derivatives 
Introducing vector notation 

D = (A, £>>.,■■•), (3.1) 
consider the KdV-type bilinear equation 

P(D„D,,Dy,--)ff^O, (3.2) 

where P is a general polynomial in D,, D^, Dy, ■ ■ ■ . 
According to ( 12.4b and ( 12.51 ), it follows that 

P(-D) = p(d). (3.3) 

Now assume that the degree of every term in P is even. 
When n is even, it's true that 



^c;;5r(-5',) /(^)/(^')|,,^, 

/=o / 

n 

j]{-i)'c'„{d:-'fix)){d'j{x)) 

1=0 
n/2-l 

2 2 (-i)'c;, {d'l-'fix)) {d'jix)) + (-i)«/2c;f {df-f{x)f 



(3.4) 



i=0 

= 2(5:;/)/ + //, 

where H = //(/^, /,■■■) is a polynomial in /,,/,■■■ , WITHOUT /, wdH^D"J f-2 {d"J) f. 
This simple observation implies the following crucial formula 

P(D„ D,, Dy, ■■■)/■/ = 2 {P(d„ d,, dy, ■■■)/)/ + ^ = 2 (P{d)f^ f + H, (3.5) 

where // is a polynomial in /i ,/,,■■■ , WITHOUT /, and P{d) - Pid,, cJ.v, dy, ■ ■ ■). Some simple 
examples are listed as follows to illustrate equation ( 13.5b 

d2/./ = 2(/,,/-/2), /7=-2/2, (3.6) 
D\f ■ f = 2(/4,/ - 4/3,/, + 3/i,), H = -8/3,/, + 6/2,, (3.7) 
D,D,/-/ = 2(/,/ -//,), // = -2//,, (3.8) 

where /t, = 5*/ = k eN. Furthermore, using equation ( 13.5b . it is easy to find P(D)f ■ / is a 
quadratic polynomial in /, /,, ft, - ■ ■ , and then to obtain 

/^ -/+/2- (3-9) 

Here Pi and Po are polynomials in /,,/, • ■ ■ , WITHOUT /. Pi = 2P(d)f, is a linear polyno- 
mial; P2 = // is a quadratic polynomial. 
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3.2. Logarithmic Transformation 
Consider the derivation formula, 

if\ fxg-fgx fx fgx 



(3.10) 



Furthermore, for arbitrary independent variable f, it's easy to show that 

{fx\ Pi Pi Pm 



(3.11) 



where A: e N and P,(/ - 1, 2, ■ • ■ , A; + 1) is a homogeneous polynomial of degree ; in /v, /?,■■•, 
WITHOUT /, and 

Pi P,^,^(-\fk\fj'^. (3.12) 

are obtained by using equation (13.101 1 repeatedly. 
Furthermore, when m + n > 1 , it shows that 

d"d"' (In/) = + ■ ■ • + (-l)"'+"-'(m + n- 1)!^^ (3.13) 

Still, by using P as the denotation of a general polynomial and letting the lowest degree of the 
terms in P be equal to or greater than 1, we have 

P(d)(lnf)^^ + ^ + --- + ^, (3.14) 

where r - degree(P) and Pi = P{d)f. Here, P,(/ = 1 , 2, • ■ • , r) also denotes a homogeneous 
polynomial of degree iin f^, f, ■ ■ ■ , WITHOUT /. 

Using the logarithmic transformation u - 2 In/, the following expressions are obtained: 



/^ _ /?\ Dlf-f 
f P) P 



(3.15) 



T / fxt fxfi \ DxDtf ■ / /a 1 



2j j.i^.'L^^ ^ -—j.\Ajx -fx ^ xf f o i ^ xf ' f 



{fxxxx ^fxxxfx + 3/2, _^ 12/,,/2 6/4 



I / P P P ) 



p \ f 



(3.18) 



Jfxxxxxx ^ 120/h Z)^/-/ ,.Dif^Dif_f (b^\' 
«6. = 2(— j = — -15 — — .30^— j . (3.19) 



Clearly, the expressions above can lead us to finding the KdV-type bilinear form associated 
with logarithmic transformation in NPDE. 
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3.3. The Relationship Between Logarithmic Transformation and KdV-type Bilinear Equation 
In this paper, we consider the logarithmic transformation 



u = 2a(ln/)„ 



n = 0,1,2 



(3.20) 



where / = f(x, f, ■ ■ ■ ) and a is a nonzero constant. 
Substituting ( IT20I ) to a NPDE 

W{u, Ux, Ui, ■ ■ ■ ) - 0, 

we get 

W(f, fx, //,■■■ ; a) = W(m, Ux, U,, ■ ■ ■ )|„=2Q-(ln /),„ 

If NPDE Wif, fx, /,,■■■; a) has the KdV-type bilinear form 



(3.21) 
(3.22) 



PiD)f ■ f 

f- 



0, 



(3.23) 



where m is a nonnegative integer, using the properties of Hirota operator and logarithmic trans- 
formation, our algorithm will find the undetermined parameters n, m, P{D), and a. 

Because of ( 13.141 ). the terms in NPDE with degree k generate the homogeneous expression 
with logarithmic transformation ( 13.201 ) 



Pk Pk+i 
+ —r-r + 



■k+l 



(3.24) 



That is, with logarithmic transformation ( 13.201 ), the terms in NPDE with degree k do NOT gener- 
ate the terms jf (i - 1,2,3, - ■ ■ ,k — I). 
Consider KdV-type bilinear form 



P{D)f-f 
P 

p 

^ 1 mx ^ 



' mx 



Pif + Pi 



J I WX 



(-l)"'(m+l)!P2/r 



(_iy"(^+ 1)1^2/™ 2&'^P{d)f) 



f 



f 



1+2 



f 



+ • ■ • + 



(3.25) 



/' 



•1+2 



Therefore, if a equation has the KdV-type bilinear form, with logarithmic transformation 
( 13.20b . the term must be generated from the linear terms in NPDE (13.21b . 

Denoting the linear part in NPDE ( 13.211 ) by P;( d )u, with logarithmic transformation ( 13.201 ). 
according to (13.14b . the linear part is transformed into 

«=2«(in/)„, 2ad'lP,{d)f 
Pi{o)u > \ . 



Thus, it follows that 



P{d)^ad"xr"'P,{d). 
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(3.26) 



From the above equation, we find it important that the linear terms of the NPDE play a crucial 
role in finding the HBF. Hence, from the linear part in NPDE, P, the function of HBP, can be 
obtained. 

Let 

/7W:/ \ ^ ^^j^ . . . . (3 
then the difference should be zero, namely. 



Res^^W(f,Uf„---;a)- 



= 0. (3.28) 

^ mx 



Clearly, « is a finite enumerated integer and m is also a finite enumerated integer (See Step 5 
of algorithm in Section 5) and P{d) = ad"-'"Piid). 

With P satisfying ( 13.26b . the coefficient of j in Res^ is equal to zero. Then the homogeneous 
expression in / and its derivatives follows as: 

Res^^^ + --- + ^. (3.29) 

Here, r is a finite positive integer 

Solving the numerators of ( 13.291 ), if there exists the nonzero numerical solution a, the unde- 
termined parameters n, m, P{D) are also achieved. It can be proven that a has no more than one 
nonzero numerical solution which satisfies (13.281) with certain n and m.The proof will be given 
in our next paper 

In summary, if a NPDE has the KdV-type bilinear form associated with the logarithmic trans- 
formation (13.20b . our algorithm will surely find it. 



4. An Example Of our Algorithm 

Now, let us take the KdV equation as an example. 

Uxxx + 6mm,; + u, — 0. (4.1) 

Here, P,id) = d, + dl- 

First substitute m = 2a; In / to ( 14.11 ) and simplify the equation, the homogeneous equation in / 
and its derivatives follows as: 

2aifxxx + 12a(ln/)/, + /,) 6a/„-/, 4a _ 

f r- p 

Since the above equation has ln(/) term, we goto the next logarithmic transformation m = 
2a(ln/),. 

Likewise, substitute m = 2(i'(ln/)^ to (14.1b and simplify the equation, the homogeneous equa- 
tion in / and its derivatives follows as: 

2(y{fxxxx + f,x) 2a(- 12a/„-/, + 4/„-,/, + 3/?, + fj,) 24a/?(-a/, + /,-,) 12a/^ _ 



From the above equation, the coefficient of — is acquired: Pi - 2a{f^xxx + fix) - 

2a {dl + dxd^f - 2adxPi(d)f . Because the order of derivative is even, the Hirota bilinear 

form can be obtained: P(D)f ■ f - a{iy\ + D^D,)/ ■ f, with which divided by and subtracts 
the above equation the difference is as follows: 

l2fxM-2afx + fxx) 24af^(-afx + /„-) 12g/f 

Solving the numerators of the difference equation, it only has a — 0. So we goto the next 
logarithmic transformation u - 2a{lnf)jcx- 

Similarly, substitute u - 2a{ln f)xx to ( 14.1b and simplify the equation,the homogeneous equa- 
tion in / and its derivatives follows as: 

2a{f,xx + fxxxxx) _ 2a(-l2afxxxfxx + fxx ft + ^fxxxxfx + '^Of XXX fxx) + 2fxftx 

f p 

4/,a(18a/2. + 6ff/„,,/, - 15/2. - 10/,„./, - /,/,) UOafxxpia - 1) 48ap(a - 1) 



P P P 



= 0. 



From the above equation, the coefficient of — is obtained: f i = 2a(fxxxxx + fxx) = 
2a{(fi^ + d\d^f. Because the order of derivative is odd, we integrate Pi and get the Hirota 

bilinear form P{D)f ■ f - a{D'l + DjcD,)f ■ f, with which divided by and derived in terms of 
X one time then subtracts the above equation, the difference is obtained: 

-24ty/„-,/,-,(g-l) 24fMa - l){3fl + /,„/,) 120g/„-/,^(g - 1) 48ap(a - D _ „ 
P ^ P P ^ P 

Solving the numerators of the difference equation, a nonzero solution, a - 1, is achieved. 

In conclusion, by using the logarithmic transformation u - 2{lnf)xx, it shows that the equation 
(14.1b is equivalent to 

i(DUDxD,)f-f\ 



P 

This example illustrates our method intuitively. 



5. An Algorithm of Bilinear 

Consider the general NPDE 



W(u, Ux, Uf,- ■ ■) — 0, 



where u - u(x, t,- ■ ■), and W is a polynomial in u and its derivatives. 
The algorithm for obtaining HBF goes in the following steps: 

Step 1 Let n -Q, f - f{x, f, ■ • • ) and a is an undetermined nonzero constant. 
Step 2 Using logarithmic transformation u = 2a (In /),„, it follows that 

W - Wif, fx, /,-■■) = W(U, Ux, M„ ■ ■ ■ )l„=2.(ln/)„ 
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Step 3 IfW has In (/) term, then goto Step 6. 

Step 4 SimpHfy the equation W to get the homogeneous expressions in / and its derivatives 

w ^1 ^2 

/ /- /' 

where P, = Pi (fx, ft, ••■)(' = 1, 2, ■ • • , r) is a homogeneous polynomial of degree / in 
/,,/,,-■■, WITHOUT/. 

Step 5 Get the coefficient of — in VV' : Pi - P\{fx,ft,' ■ ■)■ Let the lowest order of the derivative of / 

with respect to in Pi is k, and the corresponding Hirota bilinear form is P{Dx, D,,- ■ ■ )f ■ f . 
Then the following equation is obtained; 

pu.,f„---)^2d':p{f,j„---\ 

where Q <m <k and the degrees of the terms in P(fx, //,■■■) are even. That is, 

• If the lowest order derivative of / is even, then m - 0, 2, 4, ■ ■ ■ , and m < k. 

• If the lowest order derivative of / is odd, then m - 1, 3, 5, • ■ ■ , and m < k. 

Step 5.1 Let m - 0, if the lowest order of the derivative of / is even, or else let m - 1. 
Step 5.2 Integrating Pi by m times in terms of x, we get P{fx,f, ■ ■ ■) = jd^'"Piifx, f, ■ ■ ■)■ 
Hence the HBF is P(Z)^-, D„ ■■■)f -f. 

~ lP(D„D„---)f-f] 



exists, then using u - 2a (In f),,^, the NPDE is equal to ( '■ ^ I - 0- 



Step 5.3 Calculate the difference, Pe^jy = W - '■ ^ , and simplify the differ- 

\ f I mx 

ence Res^ to get the homogeneous expressions in / and its derivatives 

Pi Pr 
ReSur — ^ + ■ • ■ H . 

W p 

Step 5.4 Get the equations P, = (/ = 2, 3, ■ ■ ■ , r), and solve the equations. If a nonzero a 

(PiD,„D^--)f-f] 

Therefore, the HBF is P{Dy, D,, ■ ■ ■)f ■ f, and exit of program. 
Step 5.5 Let m = m + 2. If m < k, goto Step 5.2, or else goto Step 6. 

Step 6 Let n - n + \. If n < 2 then goto Step 2, or else the NPDE has no HBF with logarithmic 
transformation, end program and exit. 

6. Applications 

Example 1. Boussinesq equation lfl6l. 171 



Utt — Uxx - 3m^j. - M4,v =0 (6. 1) 

Solution. From the equation. Bilinear within O.s outputs a - \, u - 2\n{f) and 

{D'~-D\-D\)f-f ^ 



Example 2. Sawada-Kotera equation ifislfigll 

Solution. Within 0.032s our program outputs a = -l,u = -2(ln/)tc and 



(6.2) 



= 0. 



Example 3. Kadomtsev-Petviashvili equation (12] 

(M3.V + 6UUx + U,)jc + 3S^Uyy - 

Solution. Bilinear within 0.016s outputs a = 1, m = 2(ln/) „ and 

((Dt + D,,D, + 36'-D^)f-f 



(6.3) 



0. 



Example 4. The shallow water waves equation |l21 1 

u,dx + Ml- = 



(6.4) 



Solution. By substituting u - and using the boundary condition m,| v^co - ( 16.4b can be 
converted into the differential form 



Wxt - mxt - 'iWxWxt - 3W„W, + Wxx = 0, 

Bilinear within 0.016s outputs a - \,w - 2(ln/)v and 

i(-DlD, + Dl + DxD,)f-f\ 



(6.5) 



= 0. 



Example 5. Ito equation 12211 



r 

Utt + U^xt + 6Mj:M( + 3mM„ + llUxx I M/'^-*^ = 

xj — OO 



(6.6) 



Solution. By substituting u - Wx and using the boundary condition Ut\x^-oo - ( I6.6I 1 is trans- 
formed into the differential form 



Wxff + WAxt + ^WxxWxt + ^WxWxxt + '^WixWt = 0, 

Bilinear within 0.031s outputs a - 1, w = 2(ln/) v and 

'{DlD, + Dj)f-f] 



(6.7) 



0. 



Example 6. (2+l)-dimensional breaking soliton equation 1123 1 

r 

Uf + Puxxv + 4y6MMy + AjSux I Uydx — 

—CO 
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(6.8) 



Solution. The equation can be written as (6.9) by substituting u - Wi and using the boundary 
condition MyLv^_oo = 

Wa + PWixy + 4ySw.tW.n. + Apw^^Wy = 0, (6.9) 
Bilinear within 0.016s outputs a - j,w - ^(In/)^ and 

ll> (D.D,+l3D\D,)f- f\ 
4 P 



= 0. 



Example 7. The Bidirectional SK equation ll24ll 

5 I Uttdx + 5uxxt — 15mm, - 15m.v I m,c/x - 45m^m + 15m„m.c + 15m3.vM - u^^ — (6.10) 

Solution. Similarly, we convert ( 16.10b into the differential form by substituting m = and using 
the boundary condition m,|j:^_co = 0, m„| v^_co = 

5w„ + Swxxt - 15w.fW.H - ISWxxW, - 4-5WxxW'^ + i5w3xWxx + ^5w4xWx - W(,x — 0, (6.1 1) 

fij/wear within 0.015s outputs a - -1, w - -2(ln/) v and 

({-5D^ + D^,-5DlD,)f ■f\ 



= 0. 



Example 8. (2+l)-dimensional SK equation 12511 

2 r r 

9Mr + M5.v+ 15Mva-Ma + 15M3.vM+45MtM -5 I Myyt/x — 15MMy - 15m v I MyC/jC - 5M.f.vv = (6.12) 

Solution. With the substitution u = Wx and the boundary condition My|x^_co - 0, Myyl^^-oo = 
(I6.I2I 1 acquires a differential form 

9Wxt + W6x + 15w3.i-Wxi- + l5w4xWx + 45Wxxwl - 5Wyy - 15w.i-W.i-j - ISWxxWy - 5Wixy - 0, (6. 13) 

Bilinear within 0.032s outputs a - 1, w = 2(ln f)x and 

f (D^ - 5DlDy - 5£)2 + 9DxD,)f ■ f] 



= 0. 



Example 9. (3+l)-dimensional KdV equation 113 1 

Ut + 6UxUy + Uxxy + UAxz + 60MyM; + 10m3vM- + 20My-MA.t7 = 

Solution. Bilinear within 0.016s outputs a - — (ln/) v and 

1 {DxD, + D\D: + D\D,)f ■ / ^ 

2 f 
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(6.14) 



7. The Program Code of Bilinear 



with(PDEtools) : with(DEtools) : 

## Hirota Bilinear Method 

## Bilinear Derivative / Hirota Operator 

BD:=proc(FF,DD) local f ,g,x,m,opt; 

if nargs=l then return '*'(FF[]); fi; f,g:=FF[]; x,m:=DD[]; 

opt : =args [3 . . -1] ; if m=0 then return procname(FF,opt) ; fi; 

procnsune ( [dif f (f ,x) ,g] , [x,m-l] , opt) -procname ( [f ,dif f (g,x)] , [x,m-l] , opt) ; 

end: 

'print/BD' :=proc(FF,DD) local f,g,x,m,i; f,g:=FF[]; 

f :=cat(f , ' ' ,g) ; g:=product(D[args[i] [l]]~args[i] [2] , i=2 . .nargs) ; 

if gOl then f : = " (g) * " (f ) ; fi; f; end: 

## collect(expr,f ) ; first! 

getFnumer : =proc (df ,f ,pow: :posint : =1) local i,g,fdenom; 
if typeCdf , ) then 

g: = [op(df)] ; 

f denom: =map(denom,g) ; 

for i to nops (f denom) while f denom [i] Of ~pow do od; 
if i>nops (f denom) then Iprint (f denom) ; 

error "no term(s) or numer=0 when denom=7ol" ,op(0,f ) ~pow fi; 
g:=numer(g[i] ) ; 

if not type (expand (g) ,'+' ) then Iprint(g); 

error "Expected more than 1 term about Hirota D-operator" fi; 
return g; fi; Iprint(df); 

error "expected 1st argument be type '+'."; end: 

getvarpow:=proc(df :: function) local i ,f , var ,dif ,pow; if 
op(0,df )<>dif f then Iprint(df); error "expected diff function" fi; 
f :=convert(df ,D) ; var : = [op(f )] ; dif : = [op(op( [0,0] ,f ))] ; 
pow:=[0$nops(var)] ; f : =op(op(0 , f ) ) (var [] ) ; for i to nops(var) do 

dif : =selectremove (member , dif , {i}) ; 

pow [i] : =nops (dif [1] ) ; 

dif :=dif [2] ; 

od; pow: =zip( (x ,y)-> [x,y] , var ,pow) ; pow:=remove(has,pow,{0}) ; 
[[f ,f] ,pow[]] ; end: 

#convert to Hirota Bilinear Form 

HBF:=proc(df ) local i,c,f; if type(df,'+') then 

f:=[op(df)]; return map(procname ,f) ; fi; 

if type(df,'*') then f:=[op(df)]; 

f :=selectremove(hasfun,f ,diff ) ; c:=f[2]; f:=f[l]; 

if nops(f)<>l then Iprint(df); error "need only one diff function factor 
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f:=f[]; c:='*'(c[]); f : =getvarpow(f ) ; f:=[c,f]; return f; fi; 
if op(0,df )=dif f then f : =getvarpow(df ) ; f:=[l,f]; return f; fi; 
Iprint(df); error "unexpected type."; end: 

printHBF:=proc(PL: :list) local j ,DD,f ,C,tmp,gcdC; C : =map2 (op, 1 ,PL) ; 
gcdC:=l; if nops(C)>l then tmp : = [seqCcat (_Z , i) , i=l . .nops(C) )] ; 

gcdC:=tmp *~ C; gcdC : = ' + ' (gcdC [] ) ; gcdC:=f actor (gcdC) ; 

tmp : =selectremove (has ,gcdC ,tmp) ; gcdC:=tmp[2] ; 

if gcdC=0 then gcdC:=l fi; gcdC : =gcdC*content (tmp [1] ) ; fi; 

if gcdCol then C:=C /" gcdC; fi; DD:=map2(op,2,PL) ; 
f :=op(0,DD[l] [1] [1]); 

DD:=map(z->product(D[z[i] [l]]~z[i] [2] , i=2 . .nops(z)) ,DD) ; 
DD:=zip('*' ,C,DD) ; DD: = ' + ' (DD[] ) ; gcdC * " (DD) * cat(f,' ',f); 
end: 

## print Hirota Bilinear Transform 

printHBT:=proc(uf ,u,f,i,j , PL, alpha :=1) local DD,g,C,tmp,pl; 

pl:=printHBF(PL) ; if j>0 then print (u=2*alpha*'diff ' (ln(f) ,x$j)) ; 

else print (u=2*alpha*ln(f) ) ; fi; 

if i>0 then print( 'dif f ' (pl/f '2,x$i)=0) ; else 

print(pl/f '2=0) ; fi; NULL; end: 

guessdiff order : =proc (PL : :list,x: :naine) 

local L,minorder ,maxorder ,tmp; L:=map2(op,2,PL) ; 

L : =map(z->z [2 . . -1] ,L) ; tmp : =map(z->map2(op ,2 ,z) ,L) ; 

tmp : =map (z-> ' + ' (z [] ) , tmp) ; tmp : =selectremove (type , tmp , even) ; 

minorder:=0; if nops (tmp [1] ) <nops (tmp [2] ) then minorder:=l fi; 

tmp : =map (z->select (has , z , {x}) , L) ; tmp : =map (z->map2 (op , 2 , z) , tmp) ; if 

has(tmp,{ [] }) then maxorder:=0; else tmp:=map(op,tmp) ; 

maxorder : =min(tmp [] ) ; fi; 
if type (maxorder-minorder, odd) then maxorder : =maxorder-l fi; 
[minorder, maxorder] ; end: 

guessalpha: =proc (Res ,uf ,u, f , i , j , PL, alpha) local tmp, res , pi , f lag,k; 
flag:=l; tmp:=[op(Res)] ; tmp:=map(numer,tmp) ; 
tmp:=gcd(tmp[l] ,tmp[-l] ) ; if type (tmp, ) then 
tmp:=remove(has,tmp,f ) ; fi; if tmpOO and has (tmp , {alpha}) then 
tmp:=solve(tmp/alpha~diff order (uf) .{alpha}) ; 
if tmpONULL and has (tmp , {alpha}) then Iprint(tmp); 
for k to nops ( [tmp] ) while flag=l do 

res : =collect (expand(subs (tmp [k] ,Res)) ,f,f actor) ; 
if res=0 then pl:=subs(tmp[k] ,PL) ; 
pr intHBT (uf , u , f , i , j , pi , rhs (tmp [k] ) ) ; 
flag:=0; fi; od; fi; fi; PL; end: 

Bilinear : =proc (uf ,u, f ,x, alpha) local su,h, i , j ,gl ,CB ,PL,gdo,DD ,Res ; 
if hasf un(uf , int) then error "Do not support integral function yet. 
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Please substitute int function." fi; for j from to 2 do 
Res :=1; su: =u=2*alpha*dif f (ln(f ) , [x$j] ) ; 
h:=collect(expand(dsubs(su,uf )) ,f .factor) ; 
if hasfun(h,ln) then next; fi; 

gl :=getFnumer(h,f )/2; gl : =expand(gl) ; CB:=HBF(gl); 

gdo:=guessdiff order (CB,x) ; 

for i from gdo[l] by 2 to gdo [2] do 

if i=0 then PL:=CB; else PL : =HBF(int (gl ,x$i) ) ; fi; 

DD:=add(PL[i] [l]*BD(PL[i] [2] []) ,i=l. .nops(PL)) ; 

Res:=collect(expand(diff (DD/f ~2, [x$i])-h) ,f .factor) ; 

if Res=0 then printHBT(uf ,u,f ,i,j , PL, alpha) ; break; 

elif type (alpha, name) and has(DD, alpha) then 
Res : =guessalpha(Res ,uf ,u,f , i , j , PL, alpha) ; 

fi; od; if Res=0 then break; fi; od; PL; end: 

8. Conclusions 

To sum up, an algorithm for generating the Hirota bilinear form of NPDE with logarithm 
transformation has been proposed in this paper, and the bilinear forms of a class of NPDEs are 
obtained by using package Bilinear Then we illuminate the availability of the algorithm by 
illustrating some examples. 
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